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Reaction  turbines  basically  Francis  turbines  and  propeller/Kaplan  turbines  are  suitable  for  medium  and 
low  head  hydropower  sites.  The  management  of  the  small  hydropower  plants  is  an  important  factor,  for 
achieving  higher  efficiency  of  hydro  turbines  with  time.  Turbines  show  declined  performance  after  few 
years  of  operation,  as  they  get  severely  damaged  due  to  various  reasons.  One  of  the  important  reasons  is 
erosive  wear  of  the  turbines  due  to  cavitation.  Reaction  turbines,  however  are  more  prone  to  cavitation 
especially  Francis  turbines  where  a  zone  in  the  operating  range  is  seriously  affected  by  cavitation  and 
considered  as  forbidden  zone.  Cavitation  is  a  phenomenon  which  manifests  itself  in  the  pitting  of  the 
metallic  surfaces  of  turbine  parts  because  of  the  formation  of  cavities.  In  the  present  paper,  studies 
undertaken  in  this  field  by  several  investigators  have  been  discussed  extensively.  Based  on  literature 
survey  various  aspects  related  to  cavitation  in  hydro  turbines,  different  causes  for  the  declined 
performance  and  efficiency  of  the  hydro  turbines  and  suitable  remedial  measures  suggested  by  various 
investigators  have  been  discussed. 
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1.  Introduction 

Energy  is  the  basic  requirement  for  economic  development. 
Every  sector  of  a  country’s  economy  (industry,  agricultural, 
transport,  commercial  and  domestic)  needs  input  of  energy.  The 
installed  electricity  generating  capacity  in  India  at  present  is  nearly 
128  GW  which  is  generated  66%  through  thermal,  26%  through 
hydro,  3%  through  nuclear  and  5%  through  other  renewables. 
According  to  International  Energy  Agency  (IEA),  a  threefold  rise  in 
India’s  generation  capacity  is  expected  by  2020  [1].  As  the  non¬ 
renewable  fossil  energy  sources  continues  to  deplete,  and  realizing 
the  summits  held  at  Brazil  and  Kyoto,  to  reduce  the  greenhouse 
gas  emissions,  hydropower  has  moved  towards  the  top  power 
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development  option  to  meet  the  increasing  energy  demand.  Small 
hydropower  (SHP)  has  grown  substantially  in  the  last  10  years. 
Among  all  the  renewable  energy  sources  available,  small  hydro- 
power  is  considered  as  the  most  promising  source  of  energy.  The 
hydraulic  turbine  being  the  heart  of  a  SHP  plant  transforms  the 
potential  energy  of  water  into  mechanical  energy  in  the  form  of 
rotation  of  shaft.  Hydraulic  turbines  can  be  broadly  classified  into 
two  categories  according  to  action  of  water  on  moving  blades. 
Hydro  turbines  classification  based  on  action  of  water  over  runner 
are  impulse  turbines  and  reaction  turbines.  Impulse  turbines  are 
high  head  turbines  and  operate  at  atmospheric  pressure.  At  some 
of  hydropower  sites  where  head  is  relatively  low  and  discharge  is 
of  medium  or  high  order  it  becomes  difficult  to  operate  an  impulse 
type  turbine  as  sufficient  speed  may  not  be  obtained  to  drive  the 
end  use  machine  at  such  sites.  It  is  therefore  reaction  turbines  are 
suggested  to  be  used  and  these  considered  high  speed  turbines 
under  medium  and  low  head.  These  turbines  are  closed  turbines 


P.  Kumar,  R.P.  Saini / Renewable  and  Sustainable  Energy  Reviews  14  (2010)  374-383 


375 


operate  under  variable  pressure.  Further  these  turbines  are  more 
prone  to  cavitation  and  a  zone  in  the  operating  range  is  seriously 
affected  by  cavitation  and  considered  as  forbidden  zone  especially 
in  case  of  Francis  turbine.  The  management  of  the  small  hydro- 
power  plants  for  achieving  higher  efficiency  of  hydro  turbines  with 
time  is  an  important  factor,  but  the  turbines  show  declined 
performance  after  few  years  of  operation  as  they  get  severely 
damaged  due  to  various  reasons.  One  of  the  important  reasons  is 
erosive  wear  of  the  turbines  due  to  high  content  of  abrasive 
material  during  monsoon  and  cavitation  [2].  According  to  the 
Bernoulli’s  equation,  if  the  velocity  of  flow  increases,  the  pressure 
will  fall.  In  case  of  liquid,  the  pressure  cannot  fall  below  vapor 
pressure,  the  vapor  pressure  is  defined  as  the  pressure  at  which 
liquid  is  vaporized  at  a  given  temperature,  which  depends  upon  the 
temperature  and  height  above  mean  sea  level  of  the  site.  Whenever 
the  pressure  in  any  turbine  part  drops  below  the  evaporation 
pressure,  the  liquid  boils  and  large  number  of  small  bubbles  of 
vapors  are  formed.  It  may  happen  that  a  stream  of  water  cuts  short 
of  its  path  giving  rise  to  eddies  and  vortices  which  may  contain 
voids  or  bubbles.  These  bubbles  mainly  formed  on  account  of  low 
pressure,  are  carried  by  the  stream  to  higher  pressure  zones  where 
the  vapors  condense  and  the  bubbles  suddenly  collapse,  as  the 
vapors  are  condensed  to  liquid  again.  This  results  in  the  formation 
of  a  cavity  and  the  surrounding  liquid  rushes  to  fill  it.  The  streams 
of  liquid  coming  from  all  directions  collide  at  the  center  of  cavity 
giving  rise  to  a  very  high  local  pressure  whose  magnitude  may  be 
as  high  as  7000  atm.  Formation  of  cavity  and  high  pressure  are 
repeated  many  thousand  times  a  second.  This  causes  pitting  on  the 
metallic  surface  of  runner  blades  or  draft  tube.  The  material  then 
fails  by  fatigue,  added  perhaps  by  corrosion.  Some  parts  of  turbine 
such  as  runner  blades  may  be  torn  away  completely  by  this 
process.  The  phenomenon  which  manifests  itself  in  the  pitting  of 
the  metallic  surfaces  of  turbine  parts  is  known  as  cavitation 
because  of  the  formation  of  cavities.  Due  to  sudden  collapsing  of 
the  bubbles  on  the  metallic  surface,  high  pressure  is  produced  and 
metallic  surfaces  are  subjected  to  high  local  stresses.  Thus  the 
surfaces  are  damaged  [3].  Prof.  D.  Thoma  suggested  a  dimension¬ 
less  number,  called  as  Thoma’s  cavitation  factor  a  (sigma),  which 
can  be  used  for  determining  the  region  where  cavitation  take  place 
in  reaction  turbines: 


<7  = 


where  Hb  is  the  barometric  pressure  head  in  m  of  water,  Hs  is  the 
suction  pressure  at  the  outlet  of  reaction  turbine  in  m  of  water  or 
height  of  the  turbine  runner  above  the  tail  water  surface,  H  is  the 
net  head  on  the  turbine  in  m  of  water  [4]. 

The  value  of  Thoma’s  cavitation  factor  (a)  for  a  particular  type  of 
turbine  is  calculated  from  the  Eq.  (1).  This  value  of  Thoma’s 
cavitation  factor  (a)  is  compared  with  critical  cavitation  factor  (crc) 
for  that  type  of  turbine.  If  the  value  of  a  is  greater  than  crc  cavitation 
will  not  occur  in  that  turbine. 

For  safe  operation  (cavitation-free)  of  turbine,  it  is  evident  that: 


a  >TC 


The  following  empirical  relationships  are  used  for  obtaining  the 
value  of  Tc  for  different  turbines  [5]. 

For  Francis  turbines, 

Tc  =  431  x  10“8  x  Ns2  (3) 

For  propeller  turbines, 


T  c 


0.28  + 


1 

7.5  x  (Ns/380.78)3 


where  Ns  is  the  specific  speed  of  turbine. 


Fig.  1  shows  the  variation  of  Thoma’s  cavitation  factor  (a)  with 
efficiency.  With  increase  in  suction  head  (Hs),  a  corresponding 
decrease  in  a  is  obtained  and  this  has  no  effect  on  the  turbine 
performance  as  is  seen  from  constant  efficiency  trend.  A  stage  is 
however  reached  when  any  further  increase  in  suction  head  (Hs) 
deteriorates  the  turbine  performance  and  the  turbine  efficiency 
falls.  The  critical  value  of  cavitation  coefficient,  Tc  is  determined  by 
this  point  marking  the  change  in  the  trend  of  efficiency  [5]. 

The  present  paper  discusses  the  studies  carried  out  by  various 
investigators  in  order  to  determine  the  effect  of  cavitation  and 
identify  the  gaps  for  future  studies. 

2.  Theoretical  investigations 

Kjolle  [6]  studied  the  causes  of  damages  in  hydro  turbines  and 
found  that  the  main  causes  of  damage  of  water  turbines  were  due 
to  cavitation  problems,  sand  erosion,  material  defects  and  fatigue. 
The  turbine  parts  exposed  to  cavitation  are  the  runners  and  draft 
tube  cones  for  the  Francis  and  Kaplan  turbines.  The  effect  of 
cavitation  erosion  was  found  to  be  reduced  by  improving  the 
hydraulic  design  and  production  of  components,  adopting  erosion 
resistant  materials  and  arrangement  of  the  turbines  for  operations 
within  the  good  range  of  acceptable  cavitation  conditions. 

Duncan  [7]  provided  information  about  cavitation,  cavitation 
repair,  cavitation  damage  inspection,  date  of  inspection,  cause  of 
pitting,  runner  modifications,  cavitation  pitting  locations,  methods 
of  cavitation  pitting  repair,  and  areas  of  high  stress  in  runner. 
Typical  areas  of  cavitation  pitting  were  found  as  shown  in  Fig.  2. 

Farhat  et  al.  [8]  illustrated  the  benefits  of  the  cavitation 
monitoring  in  hydraulic  turbines  using  a  vibratory  approach.  This 
technique  was  used  in  a  large  Francis  turbine  in  order  to  validate  as 
light  modification  of  its  distributor,  which  was  intended  to  reduce 
the  cavitation  aggressiveness  and  thereby  the  related  erosion. 
Cavitation-induced  vibrations  were  measured  in  fixed  parts  of  the 
turbine  prototype  and  compared  to  those  measured  in  a  similar 
and  non-rehabilitated  turbine. 

Karimi  and  Avellan  [9]  presented  a  new  cavitation  erosion 
device  producing  vortex  cavitation.  A  comparative  study  between 
various  cavitation  erosion  situations  was  carried  out  to  verify  the 
ability  of  this  vortex  cavitation  generator  to  produce  realistic 
cavitation  erosion  with  respect  to  that  observed  in  hydraulic 
machinery.  Hardened  superficial  layers  in  specimens  exposed  to 
flow  cavitation  were  found  a  thicker  than  those  in  vibratory 
cavitation,  which  leads  to  higher  erosion  rates. 

Huixuan  et  al.  [10]  carried  out  an  on-line  monitoring  system 
for  water  turbines.  Both  audible  sound  (20  Hz-20  kHz)  and 
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Fig.  2.  Francis  turbine— typical  areas  of  cavitation  pitting. 


ultrasound  (50-300  kHz)  were  continuously  monitored.  The  signal 
characteristics  such  as  standard  deviation,  noise  level  and 
frequency  compositions  were  evaluated.  The  evaluation  results 
were  stored  in  database  in  association  with  the  operating  condition 
donated  by  water  head  and  wicket  gate  opening  (or  power  output). 
According  to  the  frequency  characteristics,  the  sound  emitted  by 
cavitation  was  distinguished  from  the  others  such  as  water  flow 
sound  and  mechanical  sound.  The  cavitation  intensity  on  different 
water  head  and  power  were  traced  out.  The  degree  of  cavitation 
erosion  was  estimated  according  to  the  cavitation  intensity  at  a 
fixed  operating  condition  such  as  rated  power  and  designed  water 
head. 

Alligne  et  al.  [11]  studied  the  hydroelectric  power  plants  for 
their  ability  to  cover  variations  of  the  consumption  in  electrical 
power  networks.  In  order  to  follow  this  changing  demand, 
hydraulic  machines  were  subjected  to  off-design  operation.  In 
that  case,  the  swirling  flow  leaving  the  runner  of  a  Francis  turbine 
may  act  under  given  conditions  as  an  excitation  source  for  the 
whole  hydraulic  system.  In  high  load  operating  conditions,  vortex 
rope  behaved  as  an  internal  energy  source  which  leads  to  the  self¬ 
excitation  of  the  system.  The  aim  of  this  paper  was  to  identify  the 
influence  of  the  full  load  excitation  source  location  with  respect  to 
the  eigen  modes  shapes  on  the  system  stability. 


3.  Experimental  studies 

Escalera  et  al.  [12]  carried  out  an  experimental  investigation  in 
order  to  evaluate  the  detection  of  cavitation  in  actual  hydraulic 
turbines.  The  methodology  was  based  on  the  analysis  of  structural 
vibrations,  acoustic  emissions  and  hydrodynamic  pressures 
measured  in  the  machine.  The  results  obtained  for  the  various 
types  of  cavitation  found  in  the  selected  machines  were  presented 
and  discussed  in  detail.  Various  types  of  cavitation  in  Francis 
turbines  were  found  are  as  shown  in  Fig.  3.  In  traveling  bubble,  the 
generalized  Rayleigh-Plesset  equation  was  found  valid  approx¬ 
imation  of  the  bubble  growth  and  it  can  be  solved  to  find  the  radius 
of  the  bubble,  ftB(t);  provided  that  the  bubble  pressure,  PB(t);  and 
the  infinite  domain  pressure,  P^t);  are  known.  The  equation  is 
expressed  as: 

PB(t)  —  Poo(f)  _  n  d2RB  3  /dRB\2  4v  dRB  2 y 

p  B  dt2  2  \dt  )  +  Rb  dt  +  pRB  1  J 


Bajic  [13]  carried  out  noise  sampling,  signal  processing  and 
analysis  and  data  processing,  analysis  and  interpretation  in  vibro- 
acoustic  diagnostics  of  turbine  cavitation.  These  were  investigated 
in  a  series  of  prototype  and  model  experiments  by  several  weak 
points  of  the  practice  and  were  identified  as  shown  in  Fig.  4. 
Improvements  and  new  techniques  were  developed.  These 
techniques  enabled  extraction  of  data  on  cavitation  details  and 
early  detection  of  detrimental  effects  met  in  turbine  exploitation.  A 
brief  review  of  weak  points  of  the  practice,  developed  improve¬ 
ments,  and  new  techniques,  as  well  as  examples  of  application  was 
presented.  Noise  power,  /(P),  at  the  turbine  power  P,  used  as  an 
estimate  of  the  cavitation  intensity,  could  be  simply  modeled  by 
the  exponential  rule: 


M 

I(P)  =  (6) 

m= 1 

where  Im  is  its  relative  amplitude  in  the  fully  developed  stage. 

Zuo  and  co-workers  [14]  carried  out  a  study  related  to 
cavitation  associated  pressure  fluctuations  in  hydraulic  systems 
which  was  an  important  phenomenon  that  affected  the  design, 
operation  and  safety  of  systems  such  as  hydropower  plant, 
pumping  station,  rocket  engine  fuel  system,  water  distribution 
network,  etc.  (Fig.  5).  The  physical  model  of  cavitation  cloud  in  the 
system  was  assumed  as  an  exciter,  which  had  its  own  character¬ 
istic  frequency  (/cav)  rather  than  simply  a  lumped  capacitance  (C) 
and,  fcav  was  a  function  of  cavitation  number,  and  operating 
condition  (Qi,  ni).  In  turbo-machinery  system,  frequency  was 
expressed  as: 

fcav  fcav  (^5  ^1 5  Ql )  C^) 

Wen-quan  et  al.  [  1 5  ]  carried  out  investigations  on  turbulent  flow  in 
a  3D  blade  passage  of  a  Francis  hydro  turbine  and  flow  was 
simulated  with  the  Large  Eddy  Simulation  (LES)  to  investigate  the 
spatial  and  temporal  distributions  of  the  turbulence.  The 
computed  pressures  on  the  pressure  and  suction  sides  agreed 
with  the  measured  data  for  a  working  test  turbine  model.  In  LES, 
large-scale  structures  could  be  obtained  from  the  solution  of  the 
filtered  Navier- Stokes  equations,  in  which  the  structures  smaller 
than  the  grid  size  need  to  be  modeled  is  expressed  as: 


dui  3 _ 

~dt  +  Wj {Ui  Uj) 


dP  1  9  UdUi  1 9tijV 

dXj  dxj  [  \  Xj  dXj )  _ 

2(i)£i2]iU]i 


9tjj 

dXj 


where  P=  p  -  &>2r2/2,zl,  is  the  large-scale  relative  velocity,  p2  is 
the  pressure  divided  by  the  fluid  mass  density,  co  is  the  constant 
angular  velocity  around  the  x3-axis,  r  is  the  radial  distance  to  the 
x3-axis,  si3k  is  a  circular  replacement  tensor,  v  is  the  kinetic 
viscosity,  and  r2J-  =  U[Uj  -  UiU]  is  the  SGS  stresses. 

Escaler  et  al.  [16]  carried  out  work  to  improve  the  cavitation 
erosion  prediction  methodology  in  hydro  turbines  by  the  use  of 
onboard  vibration  measurements  taken  on  the  rotating  shaft.  It 
was  discussed  whether  the  use  of  the  vibration  from  the  shaft  as  an 
alternative  to  the  bearing  ones  should  be  an  advantage  to  infer  the 
absolute  erosive  forces  taking  place  on  the  runner  blades. 

Athavale  et  al.  [17]  investigated  a  new  “full  cavitation  model” 
developed  for  performance  predictions  of  engineering  equipment 
under  cavitating  conditions.  All  the  test  cases  with  cavitation 
showed  plausible  results  (no  negative  pressures,  and  good 
convergence  characteristics).  Computations  on  the  water  jet  pump 
for  non-condensable  gas  concentrations  showed  sizeable  changes 
in  the  pump  head  developed.  The  vapor  transport  equation 
governing  the  vapor  mass  fraction,/,  is  given  as: 


where  v  is  the  kinematic  viscosity,  is  the  surface  tension  and 
‘p’  is  the  density. 


ft(pf)  +  V(pVf)  =  V(fV/)+Re  -  Rc 
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Fig.  3.  Different  types  of  cavitation  in  Francis  turbines:  (a)  leading  edge  cavitation,  (b)  traveling  bubble  cavitation,  (c)  draft  tube  swirl  and  (d)  interblade  vortex  cavitation. 


where  p  is  the  density,  V  is  the  velocity  vector,  r  is  the  effective 
exchange  coefficient,  and  Re  and  Rc  are  the  vapor  generation  and 
condensation  rate  terms. 

Nicolet  et  al.  [18]  carried  out  an  investigation  on  scale  model 
of  high  specific  speed  Francis  turbines.  The  analysis  of  the 


resulting  pressure  fluctuation  in  the  entire  test  rig  shows 
significant  pressure  amplitude  mainly  at  2.46  fn,  which  evidences 
the  excitation  mechanism.  A  component  of  pressure  fluctuation  at 
~2.5  fn  frequency  was  identified  all  along  the  draft  tube  walls,  the 
source  of  those  pressure  fluctuations  being  located  at  the  inner 
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Fig.  4.  The  sensors  placed  on  the  20  guide  vanes  react  cavitation  in  various  locations 
around  the  spiral. 


Fig.  5.  Cavitation  resonance  from  the  UM  Venturi. 
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Fig.  6.  Cumulative  weight  loss  of  various  industrial  alloys. 


part  of  the  draft  tube  elbow.  The  analysis  of  the  pressure 
fluctuations  phases  for  the  ~2.5  fn  frequency  provides  a  way  to 
determine  the  experimental  wave  speeds  along  the  draft  tube, 
which  were  the  key  parameters  for  a  numerical  simulation  of  the 
hydro-acoustic  behavior  of  the  test  rig.  The  simulation  carried  out 
for  the  full  test  rig,  taking  into  account  piping,  circulating  pumps 
and  the  scale  turbine  model  with  the  elbow  draft  tube  showed  that 
the  ~2.5  fn  frequency  value  corresponds  to  an  eigen  frequency  of 
the  system. 

Hart  and  Whale  [19]  carried  out  work  on  improved  weld 
surfacing  alloy  which  was  developed  and  tested  to  resist 
cavitation-erosion  in  hydro  turbines.  Typical  wear  characteristics 
experienced  in  laboratory  testing  correlated  to  actual  service 
conditions.  A  metallurgical  evaluation  showed  that  a  high 
strain,  work  hardening  austenitic  stainless  steel  produces  superior 
resistance  to  cavitational  erosion.  Several  industrial  alloys  were 
evaluated  using  the  vibratory  and  high  velocity  cavitation  test,  to 
produce  a  new  alloy  development  in  weld  surfacing.  Field  testing 
showed  an  improvement  in  cavitation-erosion  resistance  of  up  to 
800%  relative  to  308  stainless  steel.  The  cumulative  weight  loss  of 
various  industrial  alloys  are  shown  in  Fig.  6. 

Jean-Francois  et  al.  [20]  presented  a  comprehensive  research 
effort  aimed  at  understanding  the  influence  of  flow  unsteadiness 
on  leading  edge  cavitation.  Pressure  measurements  were  made  on 
the  suction  side  of  an  oscillating  hydrofoil  in  the  high  speed 
cavitation  tunnel.  Leading  edge  roughness  effects  were  investi¬ 
gated  for  both  cases  of  fixed  and  oscillating  hydrofoils.  The  cavity 
that  detaches  on  the  rough  leading  edge  was  thicker  while  its 
length  was  approximately  the  same  as  for  the  smooth  configura¬ 
tion.  Pressure  fluctuations  in  the  cavity  closure  region  were  higher 
for  the  rough  configuration.  It  was  observed  that  the  cavity 
detaching  on  the  smooth  leading  edge  at  fixed  incidence  can  be 
completely  vanish  when  the  hydrofoil  oscillates  (Fig.  7). 

The  term  cavitation  was  used  to  describe  the  phenomenon  of 
liquid-to-gas  and  gas-to-liquid  phase  changes  that  occur  when  the 
local  fluid  dynamic  pressures  in  areas  of  accelerated  flow  drop 


Fig.  7.  Stable  (left)  and  unstable  (right)  leading  edge  cavities. 


below  the  vapor  pressure  of  the  local  fluid  [21].  The  liquid-to-gas 
phase  change  was  akin  to  the  boiling  of  water,  except  that  it 
occurred  at  ambient  temperatures.  The  gas-to-liquid  phase  change 
produced  extremely  high  local  pressures  as  vapor  cavities  implode 
on  themselves.  Cavitation  commonly  occurred  in  hydroelectric 
turbines,  generally  appearing  around  guide  vanes,  wicket  gates, 
the  turbine  runner,  and  in  the  draft  tube.  Usually,  cavitation  within 
the  fluid  stream  is  not  damaging  to  the  turbine.  However,  when 
implosions  occur  near  solid  boundaries  within  the  machine,  flow 
surfaces  can  be  damaged  and  eroded.  Detection  of  the  cavitation 
phenomenon  is  straight  forward.  Large  increases  in  noise, 
particularly  in  moderately  high  frequency  ranges  (15-100  kHz) 
are  characteristic  of  cavitation.  In  addition,  vibration  levels 
generally  increase.  However,  in  a  machine  condition  monitoring 
program,  the  simple  ability  to  detect  cavitation  is  not  too 
beneficial.  The  real  need  is  to  learn  when  the  cavitation  is 
damaging  parts  of  the  machine. 

Avellan  [22]  found  out  that  design,  operation  and  refurbish¬ 
ment  of  hydraulic  turbines,  pumps  or  pump-turbine  were  strongly 
related  to  cavitation  flow  phenomena,  which  may  occur  in  either 
the  rotating  runner-impeller  or  the  stationary  parts  of  the 
machine.  The  study  presented  the  cavitation  phenomena  featured 
by  fluid  machinery  including  type  of  cavity  development  related  to 
the  specific  speed  of  machines  in  both  pump  and  turbine  mode,  the 
influence  of  the  operating  conditions,  such  as  load,  head  and 
submergence.  Therefore,  for  each  type  of  cavitation  illustrated  by 
flow  visualization  made  at  the  EPFL  testing  facilities,  the  influence 
of  cavitation  development  on  machine  efficiency,  operation  and 
integrity  were  also  discussed  in  Figs.  8-14. 

Grekula  and  Bark  [23]  carried  out  an  experimental  study  on 
cavitation  processes  in  a  Kaplan  model  turbine  and  studied  with 
the  aim  to  identify  mechanisms  that  promote  erosive  cavitation. 
The  studies  were  carried  out  with  high-speed  filming,  video  filming 
and  visual  observations  with  stroboscopic  light.  A  periodic  pattern 
of  the  cavitating  tip  vortex  was  observed.  Cavitation  at  blade  root 
observed  as  shown  in  Fig.  15.  The  main  feature  of  the  pattern  was 
that  the  cavitating  vortex  was  bent  towards  the  blade  surface  and 
transformed  into  cloud  formations.  It  was  also  found  that  these 


Fig.  9.  Typical  erosion  at  the  wall  of  the  blade  suction  side  due  to  inlet  edge 
cavitation,  shaded  area  A  of  Fig.  8. 
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Fig.  10.  Typical  erosion  at  the  wall  of  the  blade  suction  side  due  to  inlet  edge 
cavitation,  shaded  area  B  of  Fig.  8. 


Fig.  12.  Typical  erosion  at  the  wall  of  the  runner  hub  due  to  interblades  cavitation 
vortices,  shaded  area  D  of  Fig.  8. 


cloud  formations  appears  in  bands,  with  a  periodicity  which 
corresponds  approximately  to  the  spacing  of  the  guide  vanes. 
Moreover,  the  tip  vortex  cavitation  was  characterized  by  fine- 
scaled  cloud  cavitation  at  the  points  where  it  approached  the  blade 
surface,  which  indicated  an  erosion  risk. 

Harano  et  al.  [24]  investigated  two  splitter-blade-fitted  runners 
installed  in  the  Francis  turbines.  The  splitter-blade-fitted  runner 
was  a  practical  application  achieved  through  joint  research, 
utilizing  the  latest  fluid-analysis  technology,  model  testing, 
strength  analysis,  etc.,  between  Hitachi  and  The  Kansai  Electric 
Power  Co.,  Inc.  The  relationship  between  inlet  diameter  and  outlet 
diameter  of  Francis  turbine  is  shown  in  Fig.  1 6.  Results  on  vibration 
between  turbine  metal  case  and  draft  tube  are  shown  in  Fig.  17. 

Szkodo  [25]  carried  out  an  experimental  investigation  of 
cavitation  resistance  of  Fe-Cr-Mn  coating.  The  alloy  was  used  in 
mending  of  machine  elements  subjected  to  cavitation.  Chromium 
nickel  stainless  steel  0H18N9  was  used  as  the  substrate.  The 
investigated  sample  was  exposed  to  cavitation  loading  at  the 
rotating  disk  facility.  The  investigations  were  performed  in  the 
initial  period  of  the  clad  damage.  The  plastic  deformation  on  the 
investigated  surfaces  was  defined  using  an  image  analysis.  The 
microstructure,  chemical  composition  and  phase  identification  of 
the  modified  layer  were  examined  using  scanning  electron 
microscopy  (SEM),  light  microscopy  (LM),  energy  dispersive  X- 
ray  spectroscopy  (EDX)  and  X-ray  diffractometry  (XRD),  respec¬ 
tively.  The  hardness  of  processed  layers  was  investigated  by  a 


Vicker  hardness  tester.  The  results  indicated  that  there  was 
different  susceptible  to  plastic  deformation  caused  by  cavitation 
loading  for  different  kind  of  laser  processing  (Fig.  18). 

Nicolet  et  al.  [26]  investigated  is  a  1  GW,  four  Francis  turbines 
power  plant.  Francis  turbines  feature  a  cavitating  vortex  rope  in 
the  draft  tube  resulting  from  the  swirling  flow  of  the  runner  outlet. 
For  hydraulic  system  modeling  hyperbolic  differential  equations 
are 


dh  a 2  9Q.  _ 
dt  +  gA  dx~ 

dh+±d_a+wL 

9*  gA  dt  +  2 gDA2 


(10) 


(ID 


where  ‘a’  is  the  wave  speed  (m/s),  ‘Q’  is  the  discharge  (m3/s),  ‘/T  is 
the  piezometric  head  (m),  ‘D’  is  the  pipe  diameter  (m),  lA'  is  the 
cross-section  (m2),  ‘g’  is  the  gravity  (m/s2)  and  ‘A’  is  the  friction  loss 
coefficient.  The  unsteady  pressure  field  related  to  the  precession  of 
the  vortex  rope  induces  plane  wave  propagating  in  the  entire 
hydraulic  system.  The  simulation  results  revealed  the  piping 
natural  frequencies  that  were  excited  by  the  draft  tube  pressure 
source.  In  addition,  the  transfer  function  between  the  draft  tube 
pressure  source  and  the  generator  electromagnetic  torque  pointed 
out  the  risk  of  electrical  power  swing.  However,  the  risk  can  be 
really  evaluated  only  knowing  the  pressure  excitations  and  the 
draft  tube  wave  speed.  If  the  first  one  can  be  obtained  from  scaled 
model  testing,  the  latter  has  to  be  estimated,  either  experimentally 
from  vortex  rope  photography,  or  in  the  future,  by  CFD.  The 
methodology  was  a  helpful  tool  for  predicting  the  risks  of 
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Fig.  11.  “Frosted”  area  at  the  wall  of  the  runner  blade  trailing  edge  due  to  outlet 
traveling  bubble  cavitation,  shaded  area  C  of  Fig.  8. 


Fig.  13.  Limits  of  cavitation  development  within  the  operating  range  of  a  Francis 
turbine.  (1 )  Suction  side  leading  edge  cavitation  limit;  (2)  pressure  side  leading  edge 
cavitation  limit;  (3)  interblade  cavitation  vortices  limit;  (4)  discharge  ring  swirl 
cavitation  limits. 
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Fig.  14.  Typical  eroded  areas  of  a  Kaplan  runner. 


resonance  at  the  early  stages  of  pre-design  or  as  a  help  for  on  site 
diagnostic  purposes. 


4.  Analytical  investigations 

Muntean  et  al.  [27]  carried  out  a  numerical  investigation  of  the 
cavitating  flow  in  a  Francis  turbine  runner.  First,  the  steady  non- 
cavitating  relative  flow  was  computed  in  a  runner  interblade 
channel  using  the  mixing  interface  approach.  Second,  the 
cavitation  model  was  activated.  Results  for  cavity  shape  and 
extent,  as  well  as  for  the  pressure  distribution  on  the  blade  with 
and  without  cavity  were  presented  and  discussed.  The  equation  of 
pressure  coefficient  (Cp)  was  as  given  below: 


P  ~  Pref 
P^ref 


(12) 


where  pref  is  the  reference  pressure,  p  is  the  density  and  E  is  the 
specific  energy. 

Qian  et  al.  [28]  carried  out  a  three-dimensional  unsteady 
multiphase  flow  simulation  in  the  whole  passage  of  Francis 
hydraulic  turbine.  The  pressure  pulsation  was  predicted  and 
compared  with  experimental  data  at  different  positions  in  the  draft 
tube,  in  front  of  runner,  guide  vanes  and  at  the  inlet  of  the  spiral 
case.  The  relationship  between  pressure  pulsation  in  the  whole 
passage  and  air  admission  was  analyzed.  The  computational 
results  showed  that  air  admission  from  spindle  hole  decreased  the 
pressure  difference  in  the  horizontal  section  of  draft  tube,  which  in 
turn  decreased  the  amplitude  of  low-frequency  pressure  pulsation 
in  the  draft  tube;  the  rotor-stator  interaction  between  the  air  inlet 


Fig.  15.  Cavitation  at  the  blade  root. 


(min- *,kW, m) 


Fig.  16.  Relationship  between  inlet  and  outlet  diameter  of  Francis  turbine. 


and  the  runner  increases  the  blade-frequency  pressure  pulsation  in 
front  of  the  runner. 

Bajic  [29]  worked  out  novel  technique  for  vibro-acoustical 
diagnostics  of  turbine  cavitation  and  its  use  demonstrated  on  a 
Francis  turbine.  The  technique  enabled  identification  of  different 
cavitation  mechanisms  functioning  in  a  turbine  and  delivered 
detailed  turbine  cavitation  characteristics,  for  each  of  the 
mechanisms  or  for  the  total  cavitation.  The  characteristics 
specified  the  contribution  of  every  critical  turbine  part  to  the 
cavitation  intensity.  Typical  diagnostic  results  were  obtained  as 
enabled  optimization  of  turbine  operation  with  respect  to 
cavitation  erosion,  showed  it  how  a  turbine’s  cavitation  behavior 
can  be  improved;  and  formed  the  basis  for  setting  up  a  high- 
sensitivity,  reliable  cavitation  monitoring  system. 

Yuliun  et  al.  [30]  carried  out  a  study  of  cavitating  flow  in  a 
Kaplan  turbine  having  numerical  simulation  with  a  cavitation 
model  and  a  mixture  two-phase  flow  model,  which  are  in 
commercial  code  Fluent  6.1.  It  was  performed  to  unsteady 
cavitating  turbulent  flow  within  the  entire  flow  passage  of  a 
Kaplan  turbine  in  a  hydraulic  power  plant.  Based  on  the  calculation 
results,  region  and  degree  of  cavitation  occurred  in  the  flow 
passage  of  the  turbine  running  in  a  given  working  conditions. 

Muntean  et  al.  [31]  carried  out  numerical  investigation  of  the 
3D  flow  in  the  distributor  (stay  vanes  and  guide  vanes)  of  the 
GAMM  Francis  turbine.  The  domain  corresponded  to  the  dis¬ 
tributor  (stay  vane  and  guide  vane)  interblade  channel.  There  were 
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Fig.  17.  Results  on  vibration  between  turbine  metal  case  and  draft  tube. 
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Fig.  18.  Sample  surface  after  grinding  and  polishing— (a)  after  2  min,  (b)  after  4  min,  (c)  after  5  min  and  (d)  after  7  min  of  cavitation  erosion. 


three  main  issues  addressed  in  the  paper:  first,  using  the  numerical 
methodology  presented  above,  the  distributor  flow  for  several  guide 
vane  angle  values  was  computed  as  shown  in  Fig.  19.  As  a  result,  the 
guide  vane  torque  versus  guide  vane  opening  angle  was  computed 
for  the  actual  position  of  the  guide  vane  axis;  second,  it  was 
investigated  the  flow  for  the  whole  range  of  the  guide  vane  positions, 
at  four  different  locations  of  the  guide  vane  axis.  The  optimization 
criterion  considered  was  the  minimization  of  the  mechanical 
loading  of  the  turbine  regulating  system.  This  means  that  the 
extreme  value  of  the  torque  applied  to  the  guide  vane  shaft  has  to  be 
minimized.  The  designed  parameter  to  be  optimized  is  the  location 
of  the  guide  vane  shaft  axis  along  the  guide  vane  chord  line. 

Nilsson  and  Davidson  [32]  a  parallel  multiblock  finite  volume 
CFD  code  CALC-PMB  (parallel  multiblock)  for  computations 
of  turbulent  flow  in  complex  domains  was  developed  for  an 
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Fig.  19.  Three-dimensional  computational  domain  for  the  GAMM  Francis  turbine 
distributor. 


investigation  of  the  turbulent  flow  in  Kaplan  water  turbines.  This 
work  was  focused  on  tip  clearance  losses,  which  reduced 
the  efficiency  of  a  Kaplan  water  turbine  by  about  0.5%.  The 
computational  results  from  four  different  operating  conditions, 
with  different  guide  vane  angles,  were  compared  in  the  work.  It 
was  found  that,  the  computations  capture  a  vortical  structure  close 
to  the  leading  edge  tip  clearance,  where  the  tip  clearance  flow 
interacts  with  the  shroud  boundary  layer  and  cavitational  bubbles 
are  formed.  The  tip  blade  loading  increased  when  the  specific 
speed  decreased.  Detailed  measurements  need  to  be  performed, 
for  use  as  boundary  conditions  and  validation  of  the  results. 

Gagnon  and  Deschenes  [33]  focused  on  unsteady  rotor-stator 
interaction  in  a  propeller  axial  turbine.  The  flow  behavior 
was  analyzed  at  different  rotor  and  stator  relative  locations 
with  numerical  simulations  using  a  commercial  code  and  k-e 
turbulence  model.  The  main  goal  was  to  study  unsteady  flow 
phenomena  such  as  wake,  separation,  forces  and  pressure 
fluctuations  in  the  propeller  turbine.  This  investigation  will  help 
to  design  a  series  of  flow  measurements  used  in  turn  to  improve 
future  CFD  simulations  with  realistic  velocity  profiles  as  boundary 
conditions.  Wakes  and  unsteady  rotor-stator  interactions  for  a 
different  operating  regime  of  a  propeller  turbine  was  studied.  The 
wake  behind  guide  vanes  was  dissipating  very  fast.  It  was  not 
predicted  whether  it  had  an  effect  on  the  rotor  stator  unsteady 
interactions  since  higher  harmonics  of  the  spectrum  were  not 
well  defined.  The  main  interaction  was  therefore  attributed  to  the 
pressure  field  fluctuations.  Finally,  the  upcoming  experimental 
phase  to  increase  our  knowledge  and  understanding  of  the 
propeller  turbine  flow  may  be  looked  forward. 

Susan-Resiga  et  al.  [34]  carried  out  a  numerical  analysis  of  the 
swirling  flow  downstream  through  a  Francis  turbine  runner.  It  was 
found  that  the  flow  stability  characteristics  were  changed  when 
decreasing  the  discharge.  In  Fig.  20  it  was  shown  that  the  swirling 
flow  in  the  survey  section  downstream  the  runner,  in  the  draft  tube 
cone,  reached  a  critical  state  in  the  neighborhood  of  the  best 
efficiency  operating  point.  For  larger  discharge,  the  swirling  flow 
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Fig.  20.  Jet  control  technique  for  swirling  flow  in  the  discharge  cone  of  Francis 
turbines. 


was  supercritical,  and  thus  it  was  not  able  to  sustain  ax 
symmetrical  perturbations.  However,  at  partial  discharge  the  flow 
became  subcritical  and  it  was  able  to  sustain  ax  symmetric 
perturbations.  Further  investigations  revealed  that  the  axial 
velocity  and  specific  energy  deficit  in  the  central  region  were 
responsible  for  the  helical  vortex  breakdown,  also  known  as 
“precessing  vortex  rope”. 

Aschenbrenner  et  al.  [35]  used  the  mixture  model  and  full 
cavitation  model  to  compute  the  unsteady  turbulent  flow  and 
compared  with  experiment.  From  the  vaporization  and  condensa¬ 
tion  mechanics  between  vapor  and  liquid  and  the  assumption  of 
unchanged  pressure  in  cavity,  the  mass  change  rate  of  cavity  was 
proportional  to  liquid  pressure  and  vapor  pressure  difference. 

Nennemann  et  al.  [36]  studied  blade  surface  and  discharge  ring 
cavitation  with  an  emphasis  on  the  latter  using  CFD.  Localized  high 
flow  gradients  (cavitation)  as  well  as  more  global  unsteady  effects 
(rotor  stator  interaction)  contributed  to  the  phenomena  requiring 
advanced  and  resource  intensive  CFD  approaches,  notably  locally 
refined  meshes  coupled  with  unsteady  computations  including 
runner  rotation.  CFD  results  enabled  to  explain  the  presence  of 
discharge  ring  cavitation  at  distinct  circumferential  locations 
corresponding  to  the  number  of  guide  vanes  used  unsteady 


Fig.  21.  Tip  gap  and  tip  gap  vortex  cavitation  on  a  model  Kaplan  runner. 
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Fig.  22.  Cavitation  details  in  the  runner  domain  in  high  cavitating  operation. 


pressure  measurements  and  high-speed  camera  images  to 
qualitatively  support  the  validity  of  our  CFD  model.  Ultimately 
it  was  found  that  the  CFD  results  enabled  to  understand  the  physics 
underlying  the  problem  of  discharge  ring  cavitation.  Extensive 
parameters  were  studied  to  investigate  all  the  significant  influence 
factors  as  shown  in  Fig.  21. 

Balint  [37]  carried  out  a  numerical  investigation  by  computing 
the  3D  turbulent  single  phase  flow  in  the  Kaplan  turbine  runner, 
the  secondary  phase  for  the  water  (vapor)  was  activated  starting 
far  from  the  cavitation  occurrence  conditions.  The  output  pressure 
in  the  runner  domain  was  decreased  to  initiate  the  cavitating 
conditions  on  the  runner  blades.  The  unsteady  solver  with  reduced 
time  step  size  was  employed  to  capture  the  self-induced  cavitation 
in  the  runner  domain  and  3D  effects  together  with  wall-friction 
modification  onto  the  runner  blades.  The  shock  of  the  unsteady 
cavitating  flow  leads  to  high  distortion  of  the  flow  ingested 
downstream  by  the  draft  tube  of  the  turbine.  It  was  concluded  that 
unsteady  effects  of  the  flow  have  been  made  mainly  by  the 
unsteady  detachments  of  the  cavitation  at  the  blade  suction  side 
close  to  the  trailing  edge.  An  original  method  of  computing  the 
cavitation  inception  and  of  comparing  different  operating  regimes 
in  similar  cavitating  conditions  was  presented  by  plotting  the 
cavitation  volume  growth  as  shown  in  Fig.  22. 

5.  Case  studies 

Bauer  et  al.  [38]  explored  the  techniques  for  the  purpose  of 
visualizing  isolated  flow  structures  in  time-dependent  data. 
Primary  industrial  application  was  the  visualization  of  the  vortex 


Fig.  23.  Cavitation  bubbles  near  Kaplan  runner  blades. 
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Fig.  24.  Vortex  rope  in  Francis  draft  tube. 


rope,  a  rotating  helical  structure  which  builds  up  in  the  draft  tube 
of  a  water  turbine.  The  vortex  rope  can  be  characterized  by  high 
values  of  normalized  helicity,  which  is  a  scalar  field  derived  from 
the  given  CFD  velocity  data.  In  two  related  applications,  the  goal 
was  to  visualize  the  cavitation  regions  near  the  runner  blades  of  a 
Kaplan  turbine  and  a  water  pump  as  shown  in  Fig.  23.  Again,  the 
flow  structure  of  interest  can  be  defined  by  a  scalar  field,  namely  by 
low  pressure  values.  A  particle  seeding  scheme  based  on  quasi¬ 
random  numbers  was  proposed,  which  minimized  visual  artifacts 
such  as  clusters  or  patterns.  By  constraining  the  visualization  to  a 
region  of  interest,  problems  were  reduced  and  storage  efficiency 
was  gained.  They  applied  method  to  various  data  sets  from  our 
industry  partners,  visualizing  the  vortex  rope  in  the  draft  tube  of  a 
Francis  turbine  shown  in  Fig.  24,  and  cavitation  on  the  suction  side 
of  various  turbine  and  pump  runner  blades.  The  concept  of 
selectively  visualizing  relevant  flow  structures  was  proven  to  give 
additional  insight  into  their  complicated  dynamic  behavior. 

6.  Conclusions 

Cavitation  is  a  phenomenon  of  formation  of  vapor  bubbles  in 
low  pressure  regions  and  collapse  in  high  pressure  regions,  high 
pressure  is  produced  and  metallic  surfaces  are  subjected  to  high 
local  stresses.  It  is  difficult  to  avoid  cavitation  in  hydro  turbines 
which  cannot  be  avoided  completely,  but  can  be  reduced  to  an 
economically  acceptable  level.  Many  investigators  have  studied 
the  process  of  cavitation  in  hydro  turbines  through  experimental 
and  analytical  studies.  Some  of  the  investigators  have  reported  that 
in  spite  of  design  changes  in  the  turbine  components  and  providing 
different  materials  and  coatings  to  the  turbine  blades,  however  the 
improvement  in  most  cases  is  not  quite  significant.  It  is  therefore; 
further  experimental  and  theoretical  studies  are  required  for 
studying  the  impact  of  cavitation  in  hydro  turbine  to  determine 
impact  of  cavitation  at  different  values  of  parameters  which  was 
relate  to  the  cavitation  in  hydro  turbine.  CFD  based  analysis  of 
cavitation  in  reaction  turbines  could  be  cost  effective  solution  for 
an  extensive  analysis. 
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